Calculation of magnetic anisotropy energy in SmCo 5 
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SmCos is an important hard magnetic material, due to its large magnetic anisotropy energy 
(MAE). We have studied the magnetic properties of SmCos using density functional theory (DFT) 
calculations where the Sm /-bands, which are difficult to include in DFT calculations, have been 
treated within the LDA+U formalism. The large MAE comes mostly from the Sm /-shell anisotropy, 
stemming from an interplay between the crystal field and the spin-orbit coupling. We found that 
both are of similar strengths, unlike some other Sm compounds, leading to a partial quenching of 
the orbital moment (/-states cannot be described as either pure lattice harmonics or pure complex 
harmonics), an optimal situation for enhanced MAE. A smaller portion of the MAE can be associated 
with the Co-d band anisotropy, related to the peak in the density of states at the Fermi energy. 
Our result for the MAE of SmCos, 21.6 meV/f.u., agrees reasonably with the experimental value of 
13-16 meV/f.u., and the calculated magnetic moment (including the orbital component) of 9.4/xb 
agrees with the experimental value of 8.9/is. 
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The permanent magnet intermetallic compound 
SmCos has been studied extensively experimentally 1 ' 2 
and theoretically 3-11 . The interest in these materials is 
fueled by their large magnetic anisotropy energy (MAE), 
which is defined as the difference between the ground- 
state energies due to rotation of the magnetic field (mag- 
netization direction). It is generally understood that the 
main source of the large MAE in SmCos and other Sm- 
Co magnets is large magnetic single-site anisotropy of the 
Sm /-shell 9,12 " 14 . In simple terms, this means the strong 
spin-orbit coupling tries to align the Sm /-shell with the 
magnetic field, causing the /-shell to rotate with the field. 
Sm atoms in a lattice interact with the crystal field, and 
the energy of this interaction depends on the orientation 
of the Sm /-shell in the lattice. This is the leading con- 
tribution to the MAE. Note that if the crystal field is too 
small, the /-shell rotates freely with the magnetic field, 
producing no MAE, while if the spin-orbit is too small, 
the orbital moment is quenched and again no MAE ap- 
pears. We will see below that in SmCos both interactions 
are comparable, producing a large MAE. 

This basic understanding has existed for long time, 
and has been the basis of several model calculations, 
where atomic calculations for Sm have been combined 
with the crystal field parameters derived from first princi- 
ples calculations 15 . However, first principles calculations 
which could provide a quantitative analysis of different 
components of MAE in SmCos are still missing, to the 
best of our knowledge. In this paper we report such cal- 
culations, using an all-electron, full-potential, relativis- 
tic linearized augmented plane wave (FLAPW) method 
with an LDA+U extension to account for Coulomb cor- 
relations in the /-shell. 

Our analysis is organized as follows: We start by look- 
ing at the "Co part" of the MAE, the part not related 
to the Sm single-site anisotropy. We do this by investi- 
gating YC05 which forms in the same crystal structure 
as SmCos but contains no /-electrons. We calculate the 



MAE for YC05 and for the hypothetical C05 compound 
(corresponding to YC05 with Y removed) to elucidate 
the role of Y, and analyze the electronic origin of the 
MAE using the force theorem. We then move to SmCo 5 , 
where we first compute the MAE as a function of the 
Coulomb repulsion parameter U. These results demon- 
strate the nature of the Sm /-states which shows strong 
competition between the crystal field and the spin-orbit 
interactions. 

For our electronic structure calculations we have used 
the self consistent FLAPW method 16 . The Local Density 
Approximation (LDA) of Perdew and Wang 17 and the 
Generalized Gradient Approximation (GGA) of Perdew, 
Burke, and Ernzerhof 18 were used for the correlation and 
exchange potentials. Calculations were performed using 
the WIEN package 19 . Local orbital extensions 20 were in- 
cluded in order to accurately treat the upper core states 
and to relax any residual linearization errors. A well con- 
verged basis consisting of approximately 300 LAP W basis 
functions in addition to the local orbitals was used with 
the Y and Sm sphere radii set to 2.115 a.u. and the Co 
sphere radii to 2.015 a.u.. The results varied only within 
a few percent for reasonable choices of atomic radii (2.0- 
3.0 a.u.). In our previous study 11 on YC05 we estab- 
lished that the plane-wave cut-off parameters RKmax 
and Gmax of 9 and 14, respectively, are suffucient for 
MAE calculations, so the same parameters were used in 
this work. Spin-orbit (SO) interaction was incorporated 
using a second variational procedure 21 , where all states 
below the cutoff energy 1.5 Ry were included. The most 
recent version of the WIEN code, WIEN2k, includes the 
so-called p!/ 2 extension 16,22 , which accounts for the finite 
character of the wave function at the nucleus for the P1/2 
state that cannot be adequately represented as a linear 
combination of a finite number of solutions of the radial 
Schrodinger equation with 1 = 1. 

The crystal structure of SmCo 5 and YC05 is that of 
CaCu5 (PQ/mmm, No. 191). The experimental values 
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of a and c/a used in the calculation are 9.452 a.u. and 
0.792 for SmCo 5 and 9.313 a.u. and 0.806 for YCo 5 . The 
Co sites are separated into two sets of incquivalent atoms, 
Co(2c) having 2-fold multiplicity and Co(3g) having 3- 
fold multiplicity. (Figure 1) Including spin-orbit coupling 
into the calculation lowers the symmetry when the field 
lies along the plane (to Pmmm, No. 47), separating 
the 3 atoms corresponding to Co(3g) into two inequiv- 
alent sites which have multiplicities of 2 and 1, respec- 
tively. To eliminate a systematic error 11 we performed 
the calculation for both magnetization directions using 
the same, highest common symmetry group (Pmmm, in 
this case). 23 

In our previous work 11 , we found the MAE of YCo 5 
to be 0.32 meV/Co, which is 0.44 meV/Co lower than 
the experimental number of 0.76 meV/Co. 24 . What is 
important for our discussion is that both experimental 
and calculated numbers are substantially larger than the 
MAE of the hep Co (0.065 meV/Co) 25 . The question 
arises whether this enhancement is due to a different 
mutual arrangement of the Co ions, a charge transfer 
between Y and Co, or the MAE associated with the Y 
ions. To answer this question, we performed calculations 
for a hypothetical C05 compound defined as YC05 with 
Y removed with the positions of all cobalt atoms un- 
changed. We found MAE of 0.28 meV/Co (compared to 
0.32 meV/Co in YC05), proving that the reason for the 
relatively large MAE in YCo 5 is a favorable arrangement 
of the Co atoms. 

We can actually pinpoint the microscopic origin of 
this large MAE. The MAE depends on subtle differences 
in the electronic structure under rotations of the exter- 
nal magnetic field. While our actual calculations used 
two self-consistent energies for the two field directions, 
it is convenient to analyze the results using the force 
theorem 26 . According to this theorem, one can start 
from the same charge density (which is converged with- 
out spin-orbit), and then apply the spin-orbit within one 
iteration for the two different magnetization directions, 
ei and e 2 . One then sums the eigenvalues for the oc- 
cupied states in both cases, the difference of these sums 
corresponding to the MAE, good to second order in the 
change in the electron charge density. 



MAE « E? cc e,(ei) — £° cc ej(e2). 



(1) 



We used the force theorem to produce the DOS plots for 
YC05 shown in Figure 2. The two plots are practically 
indistinguishable on this scale, except near -1.5 eV, which 
shows a slight variation near the top of the peak. Also no- 
tice the smaller Co (i-band minority peak which crosses 
the Fermi level on its left side. We now plot the cor- 
responding difference in the one-electron energy (Eq. 1) 
against the number of valence electrons (Figure 3). Below 
full occupation (48 electrons), there are positive and neg- 
ative variations to the running value of the MAE which 
generally average to zero. At full occupation, 48 valence 
electrons, the MAE shows a small positive contribution 



(~ 1.5 meV/f.u.), in agreement with the results found by 
taking the energy difference of two self-consistent calcula- 
tions. This corresponds to the differences under rotation 
of the magnetization direction in the Co d-band minority 
peak which is cut by the Fermi level (Figure 2). Dopant 
atoms change the DOS and, hence, the calculated MAE 
enough so that one cannot blindly trust the results of 
a rigid shift of the Fermi energy, but this plot indicates 
that small amounts of Fe doping may initially increase 
the MAE before it decreases and changes sign, as has 
been seen experimentally 27 and theoretically 28 . 

Several difficulties are encountered going from YC05 to 
SmCos. Computationally, the main difference between 
the two compounds is that Sm includes an open /-shell 
that cannot be described in the framework of the conven- 
tional LDA theory. Indeed, uncorrected DFT (whether 
LDA or GGA) calculations incorrectly pin all of the /- 
orbitals at the Fermi energy (Ep) in SmCos and other 
systems containing unfilled /-orbitals. In order to cir- 
cumvent this problem, the previous electronic structure 
calculations 9,12-14 for SmCos did not treat the /-orbitals 
as valence electrons but as unhybridized core electrons. 
However, the width of the /-bands is ~ 1 cV, much too 
large for this approximation. Since the single-site MAE 
of Sm is due to the asphericity of the /-shells 29 , in other 
words, to the interaction of this shell with the crystal 
field, understanding the exact shape and occupation of 
the /-orbitals is crucial. Therefore we applied a Hub- 
bard U correction to the /-orbitals which naturally splits 
the /-bands into lower and upper Hubbard bands with 
nearly-integer occupancies. 

While the LDA+U method seems the proper procedure 
for handling the localized Sm /-orbitals, the question re- 
mains which form of the LDA+U should be used. There 
are several prescriptions, differing mainly in the way the 
double counted energy components are subtracted off 30 . 
The two most common are referred to in the WIEN 
code 19 as SIC and AMF (the former name is mislead- 
ing from the physical point of view, so Mazin et aP° 
suggested an acronym FLL, fully localized limit). As 
discussed by Mazin et aP° , the FLL prescription is more 
appropriate for well localized orbitals, such as /-electrons 
in SmCos, so it was used in our calculations 31 . 

The LDA+U method includes two parameters, U and 
J, that have no rigorous definition in a solid. In an 
atom, one way to define U and J is via the derivatives 
of the energy of the atomic level with respect to their 
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32. 



Jf = 



d(e/t -g/J 



(2) 
(3) 



We have applied these formulas to a quasiatom residing 
in a potential well defined by the DFT crystal potential, 
which is especially easy within the LMTO technique (the 
details of this procedure will be reported elsewhere) . This 
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gives for S111C05 Uf « 5.2 eV and J/ « 0.75 eV, which, 
as the subscripts explain, will act on the Sm /-orbitals 
in the DFT calculation. Note that this approach over- 
screens the Coulomb interaction by forcing all screening 
charge into one atomic sphere, therefore such calculated 
Uf may be slightly underestimated. 

The effect of including LDA+U in the calculation can 
be seen by looking at the spin-up (Figure 4a) and spin- 
down (Figure 4b) contribution of the Sm /-orbitals to 
the density of states (DOS). With just DFT-GGA, the 
Sm /-orbitals form a narrow band pinned at E F (not 
shown). The addition of U and J shifts the unoccu- 
pied spin- minority /-band up by approximately (U - J)/2 
and splits the spin minority /-band into the lower and 
the upper Hubbard bands, separated roughly by U - J. 
Furthermore, LDA+U enhances the Hund's second rule 
coupling 33 and tries to unquench the orbital moment by 
making the orbital moment projection m a good quantum 
number. Without spin orbit there is no Hund's second 
rule coupling, and orbital moment is quenched, that is, 
each /-band includes the same amount of m — 3 charac- 
ter as of the m = -3 character etc. Spin orbit in combina- 
tion with Hubbard U favors pure m-states. This scenario 
holds, for instance for such Sm compounds as SmAl2 and 
SmZn, where the individual spin and orbital moments are 
known experimentally 34 , and the orbital moment is be- 
tween -4.5 and -5 (1b, roughly as expected from Hund's 
second rule, in agreement with LDA+U calculation with 
a sufficiently large U. 35 . In SmCos, however, the crystal 
field is sufficiently strong to prevent total unquenching, 
and as a result the /-occupation matrix is not diagonal in 
the m representation (Table I), and the orbital moment 
is reduced to -2.8 /ib (within the calculation). 

This can be also illustrated by the DOS for SmCos, 
and comparing with the /-occupation matrix (Table I). 
The lowest peak in the DOS (Figure 4a) is clearly dom- 
inated by to = 3, corresponding to the first state in the 
Table. The next peak includes two states: a pure to = 
2 state and a nearly pure to = I state. The third peak 
is also composed of two closely lying peaks, one of char- 
acter m = 0, +2, -2, and a second of character m = +1, 
-I, +3, -3. The conduction band peak also contains two 
closely lying peaks, one with mostly m = and an ad- 
mixture of to = -2, and the second of to = -1 with an 
admixture of m = -3. The two empty states are very 
close in energy and correspond, one to m = -2 with an 
admixture of to =0, and the other to m — -3 with an 
admixture of m = -1. Note that if we apply LDA+U 
to scalar relativistic calculations, the charge state of the 
/-shell remains the same, but now the bands are formed 
by the real lattice harmonics, with the two E 2u states 
forming the unoccupied upper Hubbard bands, and, cor- 
respondingly, the E\ u , B\ u , Biui an d A.iu states forming 
the occupied bands. The distance between the lowest 
(E\ u ) and the highest (A 2u ) occupied states gives us a 
gauge of the crystal field strength: ~ 0.1 Ry. (Figure 5) 

A technical problem with LDA+U calculations is that, 
unlike conventional DFT calculation, localized orbitals 



like /-states can converge to a number of metastable 
configurations 35 . There is no guarantee that the config- 
uration shown in Fig. 5a is the true ground state. How- 
ever, we used different starting configurations but were 
never able reach a self-consistent configuration with pure 
m-states. The calculated magnetic moment, which de- 
pends crucially on the orbital configuration, agrees with 
experiment. The calculated spin moment is 12.2 [ib, and 
the Sm orbital moment is -2.8 /ib, while the Co orbital 
moments are about 0.1 /is each, leading to the total mo- 
ment of 9.9 /j,b, to be compared with the experimental 
value of 8.9 /is 36 ' 37 . As mentioned, fully unquenched or- 
bital moment would be between -4.5 and -5 hb 35 , thus 
reducing the total moment to 7.2-7.7 \ib and increasing 
the disagreement with the experiment by a factor of 2-3. 

The calculated MAE of SmCos comes from two 
sources. One is the MAE of the Co sublattice, analogous 
to that in YC05. The other is the single-site anisotropy 
of the Sm /-shell. The strong spin-orbit effect on the Sm 
/-shell is necessary due to the small MAE resulting when 
the spin-orbit effects are weak compared to crystal field 
effects, as in YC05. If there were no crystal field effects 
in SmCos, the /-states would be pure m-states, so the 
direction of the orbital quantization axis would always co- 
incide with the magnetic field, by virtue of the spin-orbit 
interaction. Assuming a crystal field interaction much 
weaker than the spin-orbit, we observe that the energy 
of the /-shell with its orbital moment aligned along 001 or 
along 100 comes from the dependence of the crystal field 
energy on the orientation of the /-shell. The stronger the 
crystal field, the larger the MAE. The fact that the calcu- 
lated /-bands in SmCos are not pure m-states indicates 
that the crystal field in this compound is strong, com- 
parable with the spin-orbit interaction. In fact, our pre- 
liminary calculations for a sister compound, Sm2Coi7 38 
indicate that the /-states there are much closer to pure 
m-states than those in SmCo 5 , that is, that the crystal 
field there is weaker, in agreement with the reduced MAE 
per Sm in Sm 2 Coi7. 39 

As discussed above, the single site MAE is defined by 
a delicate balance between the crystal field, the spin- 
orbit interaction, and the Hubbard repulsion. Since the 
LDA+U method implements the principal aspects of the 
latter, we do expect to have a reasonable description of 
MAE within this method. However, the accuracy of such 
calculations is necessarily limited. A good understanding 
of this can be gained by comparing the MAE calculated 
within two different versions of the WIEN code. Using 
the WIEN2k package 19 , in WIEN2k.01 the Hubbard cor- 
rection is applied before solving the second-variational 
spin-orbit equations, while in the WIEN2k.02 version 
it is done simultaneously. We obtain an MAE of 12.6 
meV/f.u. using the former version, and 21.6 using the 
latter. Both compare favorably with the experimental 
number of 13-16 me V/ f. u. 40-44 , emphasizing, however, 
the sensitivity of the result to the treatment of the cor- 
relation effects. 

It is worth noting that without the LDA+U correction, 
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in DFT-GGA, the calculated MAE was -11.1 meV/f.u. 
(note the wrong sign: the easy axis in the plane rather 
than axial), with the spin moment of 12.9 hb and an 
orbital moment of -1.5 hb- Also, calculations from the 
Dresden group 9 , otherwise similar to ours (we agree on 
the MAE in YC05 11 ), but with /-electrons treated as 
an "open core", produced MAE of only ~ 7meV/f.u.; 
presumably, such "atomic" description of the /-electrons 
underestimate the crystal field effects. 

Since the parameter U is not very well defined, it is al- 
ways instructive to check the dependence of the results on 
U, particularly since the above-mentioned quasiatomic 
procedure tends to underestimate U. Indeed, we found 
that if U is reduced to « 4 eV the upper Hubbard bands 
are too close to the Fermi level and the MAE is sub- 
stantially overestimated (~ 40 meV/f.u.). On the other 
hand, when U was increased to 6.0 eV the MAE, as ex- 
pected, changed little. It is also worth noting that using 
LDA 17 instead of GGA 18 improves the total magnetic 
moment (8.6 hb, with the experimental value being 8.9 
lis), but worsens the MAE (26.0 meV/f.u., to be com- 
pared with 21.6 meV/f.u. in GGA for the same value 
of U , or with the experiment, 13-16 meV/f.u.). Interest- 
ingly, both in YC05 and SmCo 5 the GGA results are in 
better agreement with the experiment for the MAE than 
LDA, although in the former compound calculations un- 
derestimate MAE, and in the latter overestimate. 

To conclude, we have performed first principle calcula- 
tions of the magnetic properties of SmCo5 using a highly 
accurate LAPW code including the LDA+U formalism. 
We obtained much better agreement with experiment 
than previous methods which treated the /-orbitals as 
open core rather than valence states. Comparing the cal- 
culation for SmCo5 with YC05 and with the hypothetical 
□Co 5 we conclude that the MAE of the Co sublatticc 
comes from a favorable arrangement of the Co atoms, 
which leads to a peak in the spin-majority DOS at the 
Fermi level. The MAE of the Sm /-shell comes from the 
interplay between the spin-orbit coupling, which tends to 
align the /-shell according to the magnetization direction 
and the crystal field, that aligns it according to the crys- 
tal lattice. In SmCo5 both interactions appear to be of 
approximately the same strength, which makes SmCo 5 
such an exceptionally hard magnet even compared with 
other Sm-Co compounds. 

This work was supported by the Office of Naval Re- 
search and DARPA Grant No. 63-8250-02. 
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FIG. 2. Calculated density of states (DOS) using the force 
theorem for magnetization along the Z-axis and along the 
X-axis in YC05. The difference is only noticeable around the 
peak near -1.5 eV. 
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FIG. 1. Crystal structure of SmCos. Shown are 2 layers of 
hexagonal Co(2c) with a Kagome lattice of Co(3g) atoms ly- 
ing between. The Sm atoms lie in the middle of the hexagons 
in the Co(2c) layers. 
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FIG. 3. Using the calculated density of states (DOS) of 
YC05 (Figure 2) for the two magnetization directions, the 
product of the energy and the difference in the DOS are plot- 
ted as a function of valence electrons/unit cell. At full occu- 
pation, 48 valence electrons, a small (~ 1.5 meV/f.u.) positive 
contribution to the MAE is related to the Co d-peak near the 
Fermi energy. 
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FIG. 4. Calculated (a) spin up and (b) spin down density 
of states (DOS) for the Sm /-orbitals with GGA and LDA+U 
(SIC-LDA, GGA) in SmCo 5 including spin-orbit. The lowest 
spin-up peak corresponds tom = +3, the second to a com- 
bination of m = +2, +1. The highest valence band peak is 
composed of m = 0, +2, -2 and m = +1, -1, +3, -3. The 
conduction band peak is composed of m = -2 and m = -3, -1. 
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FIG. 5. Calculated (a) spin up and (b) spin down density 
of states (DOS) for the Sm /-orbitals with GGA and LDA+U 
(SIC-LDA, GGA) in SmCos without spin-orbit. The occupied 
bands consist of, correspondingly, the E\ u , the B\ u and i?2u, 
and the A^ u staes while the conduction bands consist of the 
two E2u states. 



TABLE I. /-occupation matrix for SmCo 5 (17 =5 eV, mag- 
netization direction is 001). The first column lists the eigen- 
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